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行内 点 算法 


摘要 :作为 互补 问题 的 推广 , 权 互补 问题 是 一 种 重要 的 优化 问题 ,可 以 建 模 一 大 类 经 济 金融 中 的 实际 均衡 问题 。 由 于 非 
人 将 线性 优化 的 内 点 算法 推 
广 到 权 互 补 问题 。 基 于 中 心路 径 的 等 价 变换 , 提出 求解 非 负 和 象限 上 一 类 线性 权 互 补 问题 的 修正 全 牛顿 步 可 行内 点 算法 。 
在 每 次 迭代 时 ,算法 无 需 进 行 线 性 搜索 。 在 适当 假设 下 ,证 明了 算法 的 可 行 性 ,得 到 了 算法 的 迭代 复杂 度 。 数 值 实验 结果 
表明 了 算法 的 有 效 性 。 
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, A full-modified-Newton step feasible interior-point algorithm for a class of 
linear weighted complementarity problems 


WU Xinyang, ZHANG Ruijie, CHI Xiaoni, WANG Boda 
(School of Mathematics and Computing Science, Guilin University of Electronic Technology, Guilin 541004, China) 


stract: As the extension of a complementarity problem, the weighted complementarity problem is an important kind of e- 


中 


Auiilibrium problem，which could be used to model a larger class of practical equilibrium problems in economy and finance. 


8 :Because of the nonzero weight vector，the weighted complementarity problem is usually more complicated than the comple- 


ntarity problem. There is little available work about the algorithms for the weighted complementarity problem. In this 
per，an interior-point algorithm is extended from linear optimization to weighted complementarity problems. Based on an 
“equivalent reformulation of central path, a full-modified-Newton step feasible interior-point algorithm is proposed for solving 
亏 tlass of linear weighted complementarity problems over the nonnegative orthant. There is no linear search at each itera- 
tion. Under appropriate assumptions, we prove the feasibility of the algorithm, and obtain the iteration complexity. The nu- 
merical results illustrate that the algorithm is effective. 
Key words: weighted complementarity problem; full-modified-Newton step; feasible interior-point algorithm; iteration com- 


plexity 


权 互 补 问题 (weighted complementarity prob- 
lem, 简称 WCP) 是 一 类 新 的 优化 问题 ,目前 关于 权 
互补 问题 的 理论 和 算法 研究 较 少 。Potra 站 最 先 提出 
权 互 补 的 概念 , 即 找 到 一 对 属于 一 个 流 形 与 一 个 锥 的 
交集 的 向 量 , 使 得 它们 在 某 个 代数 中 的 乘积 等 于 一 个 
给 定 的 权 向 量 。 当 权 向 量 为 零 向 量 时 , 权 互补 问题 退 
化 为 互补 问题 。Potra 将 Fisher 市 场 均 衡 问题 、 二 
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次 规划 与 权 中 心 问题 转化 为 单调 线性 权 互 补 问题 ,并 
提出 了 求解 单调 线性 权 互 补 问题 的 2 种 内 点 算法 。 
之 后 ,Potrai 又 证 明了 充分 线性 权 互补 问题 的 一 些 
基本 结论 ,设计 了 一 种 校正 - 预 估 内 点 算法 。 目 前 , 光 
滑 牛 顿 法 和 内 点 算法 "是 求解 线性 优化 .互补 问 
题 01 轨 和 权 互 补 问题 "等 的 有 效 算法 ,其 中 内 点 算 
法 由 于 具有 多 项 式 时 间 复 杂 度 5 而 备 受 关注 。Koji- 
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ma 等 -给 出 了 线性 互补 问题 的 原 对 偶 内 点 算法 及 当 v 宇 0 时 , 易 证 
其 复杂 度 。Roos 等 中 首次 提出 了 线性 规划 的 全 牛顿 Xs 一 0 一 /全 =eV =ev 一 U。 
步 可 行内 点 算法 。 随 后 ,Zhang 等 "基于 修正 牛顿 方 。 则 方程 组 (2) 可 化 为 
向 提出 了 新 的 全 牛顿 步 可 行内 点 算法 。Xu 中 给 出 了 | 二 b,x 宇 0， 
线性 规划 的 基于 修正 全 牛顿 方向 的 不 可 行内 点 算法 ， iA'y 二 s 二 c,s 宇 0， (4) 
并 证 明了 该 算法 在 适当 条 件 下 具有 良好 的 复杂 度 上 xs 一 中 一 上 +， 
界 。Asadi 等 ”提出 了 求解 单调 线性 权 互补 问题 的 其 中 , 必 =(1 一 0)w,0 € (0,1)。 因 而 ,通过 求解 方 
全 牛顿 步 可 行内 点 算法 。 程 组 
鉴于 此 ,基于 Xu 中 提出 的 修正 全 牛顿 步 ,提出 AAx 一 0， 
一 种 求解 非 负 象限 上 线性 权 互补 问题 CLWCP) 的 可 | ee (5) 
行内 点 算法 。 本 算法 采用 修正 全 牛顿 步 ,避免 了 线 ee 
性 搜索 ,简化 了 算法 分 析 , 且 保证 了 迭代 点 都 位 于 可 得 牛顿 搜索 方向 (Ax ,AyyAs) 。 
中 心路 径 的 宗 邻 域内 ,最 后 分 析 了 算法 的 可 行 性 、 定义 
临 敏 性 及 复杂 度 。 数 值 算 例 结果 表明 ,本 算法 是 有 co 
效 的 。 4 一 -一 ,4 一 一 。 (6) 
AU AU 
六 互补 问题 由 式 (3)、(6), 方 程 组 (5) 可 化 为 
4d .一 0， 
之 考虑 R' 上 的 LWCP :找到 一 向 量 对 人 
(yy s,s)ER "XR"XR 使 得 | Ay td =0, (7) 
J d, +d,=(1—0e—v, 
14 .7 二 = 一 co 之 0， (1) 其 中 ,五 :=AVAX,V :一 diag(v) ,和 :二 diag(x)。 由 
CS 9 方程 组 (7) 易 知 , did, 一 0 。 
其 帅 ,AE RY ,bpER",cER ,0 ERR, yo= 定义 邻近 测度 
0 时 ,WCP 退化 为 互补 问题 。 Brossp) :=600)= |e—vl, (8) 


SZ 定义 LWCP(1) 的 严格 可 行 域 为 
是 已 = € RIXR” XR |Ax— 


b,A'y++s=c,xs >0) s 
中 心路 径 及 搜索 方向 


由 于 直接 求解 LWCP(1) 比 较 困 难 , 用 xs 一 @ = 
pe 代替 xs 一 @ 二 0, 求解 方程 组 
| 二 b,x 宇 0， 


A'y 十 s 一 c,s 宇 0， (2) 
XS — 0 一 /0e。 
假定 A 是 行 满 秩 和 矩阵 , 即 RCA)= mr。 背 内 点 


条 件 (IPO) 趾 成立, 即 (x,y,s) E F", 则 对 任意 参 


数 / > 0， 方程 组 (2) 有 唯一 解 
CC) yb) SG )，, 称 之 为 LWCP(1) 的 -中 
心 。 所 有 Ap 中 心 的 集 个 


{Cx CsyG) ss C4)) |w 盖 0) 称 为 LWCP(C1D) 的 
中 心路 径 。 当 jy 一 0 时 ,中 心路 径 的 极限 点 便 是 
LWCP(1) 的 最 优 解 。 当 xs 一 @ 过 0 时 ,定义 


_ /xs—@ 
v :二 /一 一 一 。 
从 


(3) 


引 理 19 若 w 与 v 正 交 , 则 

1 2 2 
| av | = 委 玫 | 二 |， | ae | < utrl’, 
因为 d, 与 d, 正 交 ,由 引 理 1 和 式 (8) 可 得 


1 1 2 
aad | <7latal = -vl = 


| v; 一 0)* = 
4 j=1 


| < 
el i=1 


[2 0 v0) +20|e vl +ng | 
“i=1 


[le—v ll?+20 le—v +n0]= 


a 


(9) 
同 理 可 得 


2 2 
laa < w+)", 
引 理 2” 对 任意 v € R”", 有 
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1 一 Cu) 过 过 1 十 6(v) ly No min | ( — 儿 jad,| > 
1.2 修正 全 牛顿 步 可 行内 点 算法 本 
设 刀 之 0,Gr 7,s ) E ,选择 适当 参数 0， mino) 一 gldd, -> 
ET, 使 得 
Ox,s sp ) ZT 1 5 la; d. | (12) 
a i i 本 
迭代 点 ad.1l-<G 一 0 一 So))， (13) 
x :二 x 二 Ax,y' :=y 十 Ay,s :一 5 十 As (10) 则 x Cs) o> 0 


满足 内 点 条 件 , 且 
dlx',s' ;4 ) Sr 
当 | xs 一 @ | 过。 时 ,算法 迭代 终止 。 
一 算法 1 求解 LWCP 的 修正 全 牛顿 步 可 行内 点 


) 选 择 参 数 几 盖 0,s 二 0,rE(0,1),0E(0， 

初始 点 C(x?,y 8") EF", 使 得 6 (Cx,s ;1) 去 

下 尺 :一 0。 

-2) 当 | xs 一 o | =< 

D 。 
13) 求 解 方程 组 (5) ,得 搜索 方向 (Ax ,Ay ,As) 。 


e 时 ,算法 终止 ,否则 , 转 步 


YX :一 X 十 Ax,y: 一 十 Ay， 
8 :一 8 十 As， 
/ :二 (1 一 0)7。 
一 & 十 1, 转 步骤 2) 。 


2 C 算 法 分 析 


2.1 可 行 性 和 收敛 性 分 析 


引 理 3 车 did, 
则 Gxt,y' ,sy EF', 

证 明 令 a € [0,1], 记 
Xla)=xaAx,yla)=yTaAMy,s (a)=s aAs, 
则 由 式 (3)、(6) 和 方程 组 (7) 中 第 3 式 得 


=(1—0(—6()), 


. ,LD 
eb 
GT 一 OA G 一 0)7 
已 2 
De ed 
e—@) 
154 一 2 十 < sg. 
1 二 5 
(11) 
因为 加 | 1, 由 引 理 2 


由 于 x(0) 一 x 二 0,8(C0) = 二 0,， 且 xu)， 


sl(a) 与 a 呈 线 性 关系 ,对 wa € 10,1],， 有 x (a)， 
s (a) 0。 相 应 地 , x (1),s(1) 0。 证 毕 。 
1 3 

引 理 4 令 9= ， 若 6(v) 过 三 ， 则 
时 4 5 


(x',y ,Ss ) € FY 
证 明 由 式 (9) 和 引 理 3 知 , 若 


二 GD) 十 IO) <(1 一 0)G1 一 SCo) )，(14) 


则 Cs) GE 因为 6 一 - 工 ， 所 以 式 (14) 
4Vn 
可 化 为 
1 1\? 1 
了 (6(w) + 二) [i | 
(15) 


当 Co) 去 时 , 式 (15) 成 立 , 即 (x*,y',s') EE 


& 
(V1) =v gd.. 


1 一 0 
证 明 由 式 (11) 得 


十 十 
1.2_XS 0 xDD wo 
a re CL 


引 理 6 若 (x',y',s*)E8, 则 


so_ 
eg | SCu) els 


1—0 es 
证 明 由 引 理 5 知 ， 


一 U 二 


(of ) ”一 v 


因此 ， 
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攻 3 
0 ol Se [ev a 
OWE< 
本 | 
minco) 一 一 2 全 | dd. | 1 十 ,1 一 Co) A 
由 引 理 2 可 得 由 式 (9)、(x ) 可 知 ， 
EE | 
本 x 也 
minCor) /1 一 SCo) 2 dd |。 0 了 (Oo) 十 WO) 
1 一 2 So 去 
引 理 7 若 Co 六,S) ET ， 则 2 多 
1 
1 十 .| 一 6Cu) 6 X40) 十 Vab) 
le—w ?|< lemv|1 ee | ada, | 。 (16) 
王 证 明 由 引 理 5 得 因为 6 一 _1 "|e Le 
而 A4Vn Xs 
+y2 人 
加 和 = 
a EL a .| 。 
本 注入 
© Es dCv) EE (SCo) 十 Vi0) 
© eol+trisl( -se)aal< < 
四 1 一 0 xs [0) 
二 |e—9 1+ /1-6 — (6(v) + Vn0) 
SS Je-vl 十 dd 
OS SCo) 2 (5o) ) 
引 理 8 令 和 So) 一 le 一 oo 1， 车 i . 4 
Cy ) EF, 则 4wV7 
et |e 一 Co ) ?| 1 1\? 
0 ) = 1 十 上 一 so) (so) 十 二 ) 
二 1 4 
E E 二 
CE 1+/1-50%) 0 
三 1 一 0 (17) 
OO 证明 由 9co ) 的 定义 式 (8) 得 又 因 式 (17) 的 右 端 函数 关于 SCo) 单调 增加 , 故 当 
(vd+t > 1 
Soor) = 2 dy = > i 6(v) 三 二 时 ,有 
n “2 2 
[1— C+) ， Ov) (30w) ) 
2 Te—Cw')?| = 2 ) 
1 十 mitory 全 T 二 minory。 4Vn 过 
pa 2 
因此 ,由 引 理 6 得 1 十 几 一 Co) 了 (so 7) 
; le— cw ?| 让 让 ) 
Ov < A4Vn 
CD) 
中 
| XS | -- 1 ， V2 lL 
ml 2 71 一 (3 十 辽 ) 
‘(1 a 
引 理 9 令 0 一 ,车 5(v) 过 二, 则 8(v1) 一 十 二 
An 2 ov 1 和 .y 2 
1 1 2 4 
| 
2° A4Vn 
、 ee 1 (18) 
证 明 由 引 理 4 可知, 车 8(v) 之 二 二 福 , 则 


G SEE。 由 引 理 7.8 得 


当 GCu ) 去 时 ,由 式 (16)、(17)、(18) 得 6 1) 一 
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1 生 10 个 问题 进行 求解 。 

汪 


2.2 复杂 度 分 析 


引 理 10 邻 9 二 一 一 


n 


s'—0)< 2 —0nu)’, 
二 TF 
证 明 由 vi= / 壮 二 一 得 
风 
xl sw) x so)= Ye ) [vt) wv) = 


(Cp) DE 
又 由 引 理 5 得 


,车 j) 扩 三 , 风 


CS 一 OO)TICOCT 


[(1— Owe! (wv ) 2。 (19) 


0) ye!' (v')*=(1—0)ywe! 


@\ TI 
加 一 到 要 
(1 一 0)me oz 人 e ) dd, < 


(I—Onlvl +n ld,d,|. (20) 
则 模 据 式 (12)、 (21) 和 引 理 2 可 知 ， 
Ee 0 Ce 和 一 oO) 去 
FoCo) ) 十 (1 一 0)7zmCL 一 SCo) )] < 


| 理 11 令 0= 元 ， 则 算 


,Ox ,Ss sp < 
n 


os | 
次 渗 代 ,才能 得 到 LWCP(1) 的 e- 近似 解 。 
证 明 因为 6(x',s" ;1") 牵 地 ,所 以 由 引 理 10 


得 
[2(1 — On 
. 一 OO | Ze ， 只 需 


(CCx181 —w) xs 一 0) < 
因此 要 使 x*s 


(x‘s*—0) (xs:—w)Z 
[2 —Onu =L2(— 0°nu jg 
从 而 当 2Q 一 0)*nyo 三 e 时 ,用 xs 一 wm | 过 s, 故 


渤 代 次 数 的 上 办 为 | 4Vilog 2 1 


3 数值 算 例 


为 检验 算法 1 的 有 效 性 ,运用 MATLAB 
R2016b 编程 ,在 Intel® Core i5 CPU 2.3 GHz,8.0 
GiB 内 存 ,IOS 操作 系统 的 计算 机 上 进行 数值 实验 。 
随机 生成 5 个 不 同 规模 的 LWCP(1) ,对 每 种 规模 产 


任意 选取 参数 yo 二 0, 起 始点 (x ,y" ,s") 及 权 


向 量 w 宇 0, 使 得 6(x",s";p") 达 rz。 随机 生成 行 满 
秩 和 矩阵 A 所 及 2 。 今 b 一 Ax'"， c 一 AIy 0 十 s"， 即 初 


始点 满足 内 点 条 件 。 算 法 的 终止 条 件 为 | xs 一 w | 过 
e, 记 Gap 二 上 xs 一 oO |。 在 算法 1 中 取 参 数 /一 1， 


1 


e 一 105,r 一 村 ,0 一 一 -。 表 1 为 算法 工 求解 不 同 
2 4V7z 
规模 的 LWCP(1) 的 数值 结果 ,其 中 数据 为 10 次 结 
果 的 平均 值 。 
表 1 求解 不 同 规模 线性 权 互 补 问题 的 数值 结果 
m n Iter Time/s Gap 
10 20 227 0. 022 9.537X10 °° 
20 30 283 0. 068 9.923X10° 
30 50 375 0. 221 9.710X10 
50 100 546 1. 043 9.920X10 
100 200 795 4.779 9.927X10 


例 1 考虑 Ri 上 的 LWCP(1), 其 中 


4 
1 二: 0 1 
A=| b=( 二 |3|,® 二 15|。 
0 0 1 2 
2 3 
初始 点 x (922) 8 (0352), x4" 二 1,y" 二 


(0,0)7，, 使 得 8G ,5 ;1 ) 去 3， 对 该 问题 运用 算 
法 1 进行 求解 ,迭代 78 次 后 得 到 最 优 解 ， 


x” 一 (2.477,1.477,2.000)',y” = (0. 385, 
0.499)T',s* 二 (1. 615,3.385,1.500)'。 图 1 为 迭代 
过 程 中 邻近 测度 及 Gap 的 值 。 从 图 1 可 看 出 , 随 着 
斤 测度 和 Gap ii 


的 减 小 ,邻近 


2.5 


逐渐 减 小 ,并 趋 于 0。 


2.0 


0.5 


0 01 02 03 0.4 05 0.6 0.7 0.8 0.9 1.0 
Hu 


图 1 迭代 过 程 中 邻近 测度 及 Gap 的 值 


22060000141015 _ ， 15 吴 昕 阳 .fbd 
ChinaXiv 合 作 期 和 有 
222 桂林 电子 科技 大 学 学 报 2022 年 6 月 
point algorithm based on modified Newton direction 
4 结束 语 


提出 了 求解 一 类 线性 权 互 补 问题 的 一 种 修正 全 
牛顿 步 内 点 算法 。 该 算法 基于 中 心路 径 的 等 价 变换 ， 
得 到 搜索 方向 。 证 明了 算法 经 过 修正 全 牛顿 步 后 生 
成 的 迭代 点 仍 严 格 可 行 , 且 具 有 多 项 式 时 间 复 杂 度 。 
数值 算 例 表明 了 算法 的 可 行 性 。 
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